X-Ray Absorption Analysis of NGC3516: 
Appearance of Fast Components With Increased Source Flux 



Tomer Holczer ^and Ehud Behar 



ABSTRACT 



By analyzing the X-ray spectra of NGC 3516 from 2001 and from 2006 obtained 
with the HETGS spectrometer on board the Chandra observatory, we find that 
the kinematic structure of the outflow can be well represented by four outflow 
components intrinsic to NGC 3516: -350 ± 100 km s~^, -1500 ± 150 km s~^, 
-2600 ± 200 km s-\ and -4000 ± 400 km s"^ A local component at z = could 
be confused in the spectrum with intrinsic component 3. Components 1 and 2 
have a broad range of ionization manifested by absorption from 23 different charge 
states of Fe. Component 3 and 4 are more highly ionized and show absorption 
from only 9 different charge states of Fe. However, we were able to reconstruct the 
absorption measure distribution (AMD) for all four. The total column density 
of each component is A^^ = (1.8 ± 0.5) x 10^^ cm-^ (2.5 ± 0.3) x 10^^ cm'^, 
(6.9 ± 4.3) X 10^2 cm-2, and (5.4 ± 1.2) x 10^^ cm'^, respectively. The fast 
components 3 and 4 appear only in the high state of 2006 and not in 2001, while 
the slower components persist during both epochs. On the other hand, there is 
no significant absorption variability within days during 2001 or during 2006. We 
find that covering factor plays a minor role for the line absorption. 

Subject headings: galaxies: active — galaxies: individual (NGC 3516) — tech- 
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INTRODUCTION 



Approximately half of type 1 AGNs sho w complex absorption in the soft X-ray band 
( ICrenshaw et al.l l2003l : iPiconcelli et al.l |2005| ) . T he blue s h ifted absorption lines come from 
a highly or partly ionized outflow first noted by iHalperru ( ll984j ). These outfiows may play 
a central role in cosmological feedback, in the metal enrichment of the intergalactic medium 
(IGM) and in understanding black hole evolution. However, the outfiow physical properties 
such as mass, energy and momentum are still largely unknown. A necessary step to advance 
on these issues is to obtain reliable measurements of the ionization distribution and column 
densities of the outfiowing material as well as the location of these outfiows, in order to 
determine the mass outfiow. 



NGC 3516 is a Seyfert 1 galaxy at a r edshift z = 008836 flKeell 119961). It exhibits 



of a few narrow components (iTurner et al 



a variable relativistic Fe Ka e mission line ( Nandra et al.l 119971 . Il999[ ) which could consist 



2002). UV observations of NGC 3516 showed 



several outfiowing components. iKriss et al.l ( 119961 ) found a low velocity component of ~ - 
300 km s~^ while also o bserving lines at close to the systemic redshift velocity (-2649 km s~^). 
Crenshaw et al.l (119991 ) revealed four slow cora ponents with outfiowing velocities of -31, -88, 
-148, and -372 km s~^. iKraemer et al.l ( 12002| ) later detected UV absorption up to velocities 
of - 1300 km s-^ 

There are many studies in the literature of the ionized X-ray absorber of NGC 3516, the 
most important of which are listed in Table [H Not all reports on the outfiow of NGC 3516 
are consistent, as described below. iKolman et al.l ( 119931 ) analyz ed a Ginga obs e rvatio n from 
1989 and found an absorber with a partial covering factor. iMathur et al.l (119971) found 



a -50 km s ^ outflowing component using a ROSAT observation from 1992. iKriss et al. 



(119961 ) fltted ASC A data fro r n 1994 with two absorption components, both outflowing with 



< -120 km s ^. iReynoldsl (119971 ) fitted the 1995 ASCA data with a single component. 



Costantini et al.l (120001 ) studied BeppoSAX data from 1 996 and 1997 and fi tted two ionization 



components outfiowing at velocities of -500 km s~^. iNetzer et al.l (120021 ) fitted ASCA data 
from 1998 and Chandra LETG data from 2000 with a single compo nent. Using 2001 pa rtially 
simultaneous Chandra, BeppoSAX, and XMM-Newton observations, iTurner et al.l (120051 ) iden- 
tified three components; The -200 km s"^ component was called the "UV absorber", the 
second was called "High" at -1100 km s~^, and the third that was called "Heavy", was 
also at -1100 km s~^, but coverage of the so urce varied, which explained changes in the 
continuum curvature. iMarkowitz et al.l (120081 ). based on a Suzaku 2005 observation found 
two component s, one was fitted wit h an outfiow velocity of -1100 km s~^, while the other 



was stationary. iTurner et al.l (120081 ) studied Chandra and XMM-Newton observations from 
2006 and identified four components. Two components were static while the third ("heavy" 
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which also had a partial covering factor) and fo urth ("high") had higher outflow velocities of 
-1575 and -1000 km s~^, respectively. Recently, iMehdipour et al.l ( 120101 ) analyzed an XMM- 
Newton observation of NGC 3516 from 2006 and identified three components outflowing at 
velocities of -100, -900 and -1500 km s~^. The fastest component also had a partial covering 
factor. 

In this paper, we wish to provide a definitive investigation of the physical conditions 
in the NGC 3516 X-ray absorber, using the archival HETGS observation from 2001 and 
2006 and exploiting its high spectral resolution, but with special focus on the ionization 
distribution of the plasma, and on lines that do not seem to fit the simple two-velocity 
picture. While the majority of studies on the X-ray spectra of AGN outflows empl oy gradually 



increasing number of ionization components, until the fit is satisfactory (e.g., iKaspi et al. 



200ll : ISako et al.ll2003l ). it is instructive to reconstruct the actual distribu tion of the column 



densi ty in the plasma as a continuous function of ionization parameter f (ISteenbrugge et al. 



20051 ). which we termed the absorption measure distribution {AMD, iHolczer et al.l 120071 ). 
Although for high-quality spectra such as the present one, two or three ionization components 
might produce a satisfactory fit, the AMD reconstruction is the only method that reveals the 
actual distribution including its physical discontinuities (e.g., due to thermal instability), and 
ultimately provides a more precise measurement of the total column density. Furthermore, 
the AMD shape can impos e tight constraints on physical outflow models and on the d ensity 
distribution in the outflow Jchelouchell2008l : iBehaiiboOol : IPukumura et al.lboiol . boilh . 



2. DATA REDUCTION 

NGC 3516 was observed by Chandra/UETGS on 9-10 April and on 11 November, 2001, 
as well as on 9-14 October, 2006 for a total exposure time of 407 ks. In each observation, 
the first three HEG (high energy gratings) and MEG (medium energy gratings) refraction 
orders (±1, 2, 3) were reduced from the Chandra archive using the standard pipeline software 
(CIAO version 4.1.2). The total number of counts in the first order (±1) between 2 and 
25 A is 141663 for MEG and 75043 for HEG. More details on the observations can be 
found in Table. [2l No background subtraction was required as the background level was 
negligible. Flux spectra were obtained by first co-adding count spectra and convoluting 
with the broadest line spread function (MEG first order) to ensure uniformity. The total 
count spectra were then divided by the total effective area curve (summed over orders) and 
observation time. Finally, spectra were co rrected for neutral Galactic absorption of Nh = 



3.23x10^° cm-2 fiDickev fc Lockman]ll990l ). Spectra from April 2001, November 2001 and 



October 2006 are presented in Fig. [T] with a binning of 50 mA. 
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SPECTRAL MODEL 



Variations of approximately 40% on time scales of few ks were observed during both the 
2001 and 2006 observations of NGC 3516. The average X-ray flux of NGC 3516 varies greatly 
over the span of the 17 years it has been observed, as can be seen from the continuum flux 
levels quoted in Table [TJ The difference between minimum and maximum reaches a factor 
of ~ 10-30 on time scales of a few r nonths, in co n trast with the impression given by the 
sequence of observations reported in iNetzer et al.l (120021 ) of slow decay in flux over a few 
years. The present work deals with the mean properties of the ionized absorber. As a first 
step we reduce a spectrum from each separate observation. It can be seen from Fig [T] that 
on a time scale of a day, the average flux varies by a factor of a few, during both the 2001 
observations (from black to red) and the 2006 observations (from cyan to yellow). 

Since the 2001 flux is consistently lower than that of 2006, we combine the three ob- 
servations from 2001 (low state) separately from the five 2006 observations (high state). 
Combining the spectra was carried out by first adding all of the counts and only then divid- 
ing them by the time-weighted sum of the effective area and total exposure time. The flux 
level in the 2006 observation was on average five times higher than in 2001. The high and 
low state spectra are shown in Fig. |2] with a binning of 50 mA. 

It can be seen that the 2006 spectrum shows much better resolved absorption troughs 
and was used in order to obtain th e warm absorber parameters. The fitting procedure 



follows our ion-by-i on fitting method ( iBehar et al.ll200ll : ISako et al.ll2001t iBehar et al.l 12003 



Holczer et al.ll2005l ). First, we fit for the broad-band continuum. Subsequently, we fit the 
absorption features using template ionic spectra that include all of the absorption lines and 
photoelectric edges of each ion, but vary with the broadening (so-called turbulent) velocity 
and the ionic column density. Covering factor of unity is used throughout this process. 
The "black" troughs of the leading lines of O"*"^ and O"^'' strongly support this assumption. 
Strong emission lines are fitted as well. The emission lines are added after the absorption 
components were modelled (hence, the emission lines are not absorbed). Another interesting 
feature that can be seen in Fig [2] is the softening of the s pectrum at higher flux levels. 
This was already observed for NGC 3783 ( INetzer et al.l 120031) and in fact expe cted from the 
cooling of a comptonizing corona above the accretion disc (iHaardt et al.ll200l[ ). 



3.1. Continuum Parameters 



The continuum X-ray spectrum of most AGNs can be characterized by a high-energy 
power-law and a soft excess that rises above the power-law below ~1 keV. This soft excess is 
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often modeled with a blackbody, or modified blackbody, altliougli it clearly is more spectrally 
complex and possibly includes atomic features. For the 2006 spectra, we used a power law 
with a photon spectral index of F = 1.48, which was fitted to the 2-6 A band, and a blackbody 
temperature of kT = 110 eV. This is a rather flat slope, which could be due to the band in 
which it is fitted, but it still provides an good fit to the spectrum (see Fig. H]). 



3.2. The Ionized Absorber 

The intensity spectrum /ij(z/) around an atomic absorption line i j can be expressed 



as: 



/o(z/) e" 



(1) 



where Io{i^) represents the unabsorbed continuum intensity, Uijiv) denotes the line absorption 
cross section for photo-excitation (in cm^) from ground level i to excited level j. If all ions 
are essentially in the ground level, Nion is the total ionic column density towards the source 
(in cm~^). The photo-excitation cross section is given by: 



Tce 



nipC 



(2) 



where the first term is a constant that includes the electron charge e, its mass mg, and the 
speed of light c. The absorption oscillator strength is denoted by fij, and 0(z/) represents 
the Voigt profile due to the convolution of natural (Lorentzian) and Doppler (Gaussian) 
line broadening. The Doppler broadening consists of thermal and turbulent motion, but in 
AGN outflows, the turbulent broadening is believed to dominate the temperature broaden- 
ing. The Natural broadening becomes important when the lines saturate as in our current 
spectrum, e.g., the O"*"® line. Transition wavelengths, natural widths and oscillator strengths 
were calculated using t he Hebrew University Lawrence Livermore Atomic Code (HULLAC, 
Bar-Shalom et al.ll2001 ). Pa r ticularly important for AGN outflows are the inner-shell absorp- 
tion lines ( Behar et al.l 2001 : Behar fc Netzer 2002). Mor e recen t and improved atomic data 



Gu et a. 



(120061) . These atomi c wavelengths 



for the Fe M-shell ions were incorporated from 
were tested against t he HETG spectra of NGC3783 fiBehar fc Netzerl l2002l : iNetzer et al. 



20031 : IGu et al.l l2006l : iHolczer et al.l 120071 ) and found to agree with the data to within the 



instrumental precision. 



Since the absorbing gas is outflowing, the absorption lines are slightly blue-shifted with 
respect to the AGN rest frame. Although blue shifts of individual lines can differ to a small 
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degree, we can identify two kinematic components with best-fit outflow velocities of f = 
—350 ±100 km s~^ and v = —1500 ±150 km s~^ which we call component 1 and component 
2, respectively. We also identify two high ionization components a.t v = —2600 ±200 km s~^ 
and V = —4000 ±400 km s~^ (component 3 and component 4, respectively). The velocities 
are set in the model to one value (for each component) for all of the ions. Figure. [3] show 
the absorption troughs of Si and Mg K-shell ions as well as Fe"*"^^, Fe+^^, Fe"*"^^ and Fe^^^ in 
velocity space, where the vertical blue lines represent the best-fit velocity components. The 
first two components are detected in all ions while component 3 has more subtle clues, like 
a "knee" in the troughs of Si^^^, Si^^^ and Mg"^^^, or a shallow trough in Mg+^'' and Fe"*"^^, 
and a much clearer trough in Fe"*"^^. Component 4 is even harder to detect. There are signs 
for it either as a small "knee" or very small troughs. The broad Ka troughs of Fe+^^ and 
Fe"*"^^, not shown in figure [3l also extend from ~ 0-4000 km s~^. However, this result can not 
be too meaningful, as the instrumental resolution at these lines is worse than 3000 km s~^. 
Fe"*"^^ does not show absorption from the faster components 3 and 4, and we conclude these 
two components appear only in higher ionization states. Above ~ 20 A there are absorption 
lines mostly from oxygen which appear to have a high velocity component matching the 
third component (—2600 km s~^), but could also be local (z=0), which is discussed in detail 
in m 

For all of the four components, a Doppler turbulent velocity fturb = 300 km s^^, referred 
to by some as b[= v^a =FWHM/(2\/h^)] is used (i.e., full width half maximum FWHM = 
500 km s^^), which is approximately the MEG broadening (23 mA) at 14A. This value 
provides a good fit to the strongest absorption lines in the spectrum. Finally, the model 
includes also the 23 mA FWHM instrumental broadening (convolved w i th the model). The 



value of Wturb = 300 km s ^ is the same value used by [Turner et al.l (120051 ) for the slow 



compo nent. Much narrower lines can not be resolved in the pr esent spectrum. [Turner et al. 



(120081 ) used Vmrh = 200 km s ^ for the slow components while iMehdipour et al.l ( l2010l ) used 



i^turb = 50 km s ^ for the slowest component and fturb =400 km s ^ for the second one. 

Our model includes all of the important lines of all ion species that can absorb in 
the HETGS waveband. In the different components of NGC 3516, we find evidence for the 
following ions: Fe+^- Fe+^s, all oxygen charge states, Ne+^- Ne+^ Mg+'^- Mg+'^\ Si+^- 
Ar+16- At+", Ca+i«- Ca+^^, and S+^^- S+^^ We also include the K-shell photoelectric 
edges for all these ions although their effect here is largely negligible. When fitting the data, 
each ionic column density is treated as a free parameter. A preliminary spectral model is 
obtained using a Monte-Carlo fit applied to the entire spectrum. Subsequently, the final fit 
is obtained for individual ionic column densities in a more controlled manner, which ensures 
that the fit of the leading lines is not compromised. Ionic column density uncertainties are 
calculated by varying each column density (while the other ions are fixed) until Ax^ = 1, as 
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described in more detail in iHolczer et al.l (120071 ). 

The best fit model is plotted over the data in Fig. |H It can be seen that most ions 
are reproduced fairly well by the model. Note that some lines could be saturated, e.g., the 
leading lines of O"*"^ and O'''''. In these cases, the higher order lines with lower oscillator 
strengths are crucial for obtaining reliable Nion values. 



3.3. AMD Method 

The large range of ionization states present in the absorber implies that the absorption 
arises from gas that is distributed over a wide range of ionization parameter ^. Throughout 
this work, we use the following convention for the ionization parameter ^ = L/(n//r^) in 
units of erg s~^ cm, where L is the ionizing luminosity, nn is the H number density, and 
r is the distance from the ionizing source. We apply the Absorption Measure Distribution 
{AMD) analysis in order to obtain the total hydrogen column density along the line of 
sight. The AMD can be expressed as: 

AMD = dNH/ d{\ogi) (3) 

and 

Nh = j AMD rf(logO (4) 
The relation between the ionic column densities Nion and the AMD is then expressed as: 

N^on = aJ ^^/_(logOrf(logO (5) 

where Nion is the measured ion column density, A^ is the element abundance with respect to 
hydrogen assumed to be constant throughout the absorber, and /jon (log is the fractional 
ion abundance with respect to the total abundance of its element. We aim at recovering the 
AMD for the different kinematic components of NGC 3516. 

For the AMD, we seek a distribution dNu / d{\og^) that after integration (eq. E]) will 
produce all of the measured ionic column densities. In this procedure , one must take into 



accou nt the full dependence of fion on ^. We employ the XSTAR code (IKallman &: Bautista 



20011 ) version 2.1kn3 to calculate /ion (log using the continuum derived in §3.H and ex- 



trapolated to the range of 1 - 1000 Rydberg. During the fit, the AMD bin values are the 
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only parameters left free to vary. The AMD errors are calculated by varying each bin from 
its best-fit value while the whole distribution is refitted. This procedure is repeated until 
Ax^ = 1. The fact that changes in the AMD in one bin can be compensated by varying 
the AMD in other bins dominates the AMD uncertainties. This is what limits the number 
of bins and the AMD resolution in ^ or T. Indeed, we choose the narrowest bins (in log^) 
that still give meaningful errors. AMD in neighboring, excessively narrow bins can not be 
distinguished by the data, i.e., different narrow-bin distributions produce the measured Nion 
values to within the errors. More det ails on the AMD binning method and error calculations 
can be found in iHolczer et al.l ( 120071 ) . 



The current method obtains a well defined distribution of ionization, which is tightly 
constrained by the data, instead of the more traditional method, which constructs a superpo- 
sition of individual ionization components that are 6 functions of C,, and are obtained through 
a global-fitting procedure. The current method should not be viewed as an upgraded version 
of the traditional one. In fact, the two approaches are fundamentally different. While the 
AMD is a bottom-up approach that uses directly measured quantities Nion to derive the 
distribution, global fitting uses a top-bottom method that imposes a physical model and 
obtains its best-fit parameters. 

The advantage of the AMD method is that it helps i dentify column- density and ioniza- 
tion trends that can then be compared with models (e.g., iFukumura et al.l 120101 ). as well as 



the precise temperature boundaries between different phases of the possi bly multi-phase ab- 
sorber, e.g., due to thermally unstable temperatures (IHolczer et al.ll2007l ). It is not intended 
to, and indeed does not necessarily provide a superior statistical fit to the spectra. 



3.4. Narrow Emission Lines 

The present NGC 3516 spectrum has a few narrow, bright emission lines. The emission 
component is added only after the continuum is set and the absorption fit is completed. 
0+^ Ka and Ne+^ Ka forbidden lines are well fitted with a simple Gaussian broadened by 
235 km s-^ FWHM (a = 100 km s"^) while Fe Ka lines were fitted with a broader Gaussian 
of 3500 km s"^ FWHM (cr = 1500 km s"^). All lines are found to be stationary to within 
^100 km s~^. The centroid wavelength and photon flux are measured for each line and 
listed in Table O 
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4. RESULTS 

4.1. Ionic Column Densities 

The best-fit ionic column densities are listed in Table [Hand the resulting model is plotted 
over the data in Fig. [H The Errors for the ionic column densities were calculated in the same 



manner as m 



Holczer et al.l ( 120071 ) . It can be seen that both component 1 and component 
2 have a wide range of ionization parameter. They both have similar absorption with ionic 
column densities of the order of 10^^ to 10^'' cm~^ in iron, silicon, neon, and magnesium, 
while those of the more abundant O ions are higher and reach ~ 10^^ cm~^. Components 
1 and 2 are similar in the amount of absorption, even though component 2 seems to have 
slightly higher column densities. Component 3 and component 4 consist exclusively of high 
ionization species, mostly K-shell ions of Ne, Mg, Si, Fe, S, Ar and even Ca. In Fe, we detect 
only the highest ionization L-shell ions in these fast components (See Figure [3D • 



4.2. AMD For Components 1 and 2 



The best-fit AMD for components 1 and 2 in NGC 3516 is presented in Fig. [5] and the 
integrated column density is presented in the bottom panel of Fig. |5l These AMDs were 
obtained using 21 charge states of Fe from Fe+^ through Fe+^^. K-shell Fe is heavily blended 
for all the kinematic components (1 through 4), however most of the absorption is due to 
components 3 and 4. Many M-shell ions are only tentatively detected. Nonetheless, the 
strict upper limits on so many Fe ions provides tight constraints on the AMD around where 
these ions form. The result of Fig. 5 is enabled by the AMD method and, obviously, can 
not be as well quantified with a standard multi - ^ fit. The current integrated AMD of the 
absorber in NGC 3516 (Fig. [5]) gives a total column density of Nh = (1.8 ± 0.5) x 10^^ cm~^ 
for component 1 and Nh = (2.5 ± 0.3) x 10^^ cm~^ for component 2. 

Both AMDs feature a statistically significant minimum at 0.7 < log^ < 1.5 (erg s^^ cm), 
which corresponds to temperatures 4.5 < logT < 5 (K). This discontinuity in the AMD 
at th e same temp eratures was also ob served in IRAS 13349 +2438. NGC 3783 (iHolczer et al. 
2007h . NGC 7469 felustin et al1l2007h . and MCG -6-30-15 faolczer et aUboioh. Indee d, the 



multi-phase nature o f AGN outflows is found in man y warm absorbers (ISako et al.l 12001 



Krongold et al.ll2003l . |2005| . 120071 : iDetmers et al.ll201l[ ). It is mostly a manifestation of the 
relatively low ionic column densities observed for the ions Fe^^^- Fe"*"^^, as can be seen in 
Table HI One way to explain this gap is th at this temperature regime is thermally unstable 
( jHolczer et al.l 120071 : iGoncalves et al.ll2010l ). Gas at 4.5 < logT < 5 (K) could be unstable 
as the cooling function A(T) generally decreases with temperature in this regime. Such 



- 10 - 



instabilities c ould result in a mu lti phase (hot and cold) plasma in pressure equilibrium, as 



suggested by iKrolik et al.l ( Il98ll ) . 



It is somewhat surprising that the AMDs of component 1 and 2 are so similar, which 
implies the physical conditions in component 1 and component 2, the distance from the 
source and the outflow density, as well as the overall outflow column are similar. Given the 
similar AMDs we are led to think that these components are connected. What about the 
faster components 3 and 4? Are they also connected to component 1 and 2?. Components 
3 and 4 have yet to be reported, so no comparisons can be made with other authors. These 
new components and a plausible geometric explanation are addressed in detail in §4.7[ while 
their ionization distributions are discussed below. 



4.3. AMD for components 3 and 4 

The best-fit AMD for component 3 and 4 in NGC 3516 is presented in Fig. O and 
the integrated column density is presented in the bottom panel of Fig. [61 The black line 
represents component 3 while the blue dotted line represents component 4. These AMDs 
were obtained using L-Shell and K-Shell charge states of Fe only as no M-shell Fe ions are 
detected, which is why these AMDs begin at log^ = 2 (erg s~^ cm). On the other hand, the 
AMD value of these component for log.^ > 2 (erg s~^ cm) is an order of magnitude higher 
than in components 1 and 2. In the K-shell Fe lines, all the components are blended, however 
the troughs are dominated by components 3 and 4. It can be seen that both AMDs show 
similar shape but component 4 has a minimum at 2.8 < log.^ < 3.4 (erg s~^ cm), mostly 
due to the low columns of Fe+^°-Fe"'"^^, as can be seen in Table |H Similarly to components 
1 and 2, we are led to think that component 3 and 4 are also connected; Particularly, since 
they both appear in the 2006 spectrum, but are absent in the 2001 spectrum, together (see 
TO]) . The total column density (integrated AMD) gives A^^ = (6.9 ± 4.3) x 10^^ cm'^ 
for component 3 and A^^^ = (5.4 ± 1.2) x 10^^ cm~^ for component 4. 



4.4. Total Column Density 

In order to further compare our results with previous outflow models for NGC 3516, we 
can formally rebin the AMD in Fig. |5]to two regions, one below log^ < 0.5 and one above 
log(^ > 1.5 (this process is done to each kinematic component). The physical parameters of 
these two ionization regions are subsequently compared with all previous works in Table O 

The total column density that we find for the low-velocity -350 km s~^ component 1 is 
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roughly 2 x 10 cm and is comparable to that of lKriss et al. Jl996h (their co mponents 1 and 



2). Most other works are consistent with this result to with i n a factor of 2 (IRevnolds 



Mathur et all 119971: iNetzer et al. 120021 : iTurner et al.l l2005l . 120081 : iMehdipour et al 



1997 



2010|). 



However, if the ICostantini et al.l ( 120001 ) hot component also refers to our component 1, then 
there is a larger disagreement of around order of magnitude. 

Component 2 (v=-1500km s~^), was detected only with the more recent grating in- 
struments. The co lumn density we find for c omponent 2 is in good agreement with the 
Turner et al.l (120051 ) " High" component and the iMehdipour et al.l (120101) component " B" and 



still consistent with the weak constraints of Markowitz et al. 



(l2008h . It is an order of mag - 



nitude lower than iTurner et al.l ( 120081 ) "Zone 3" and "Zone 4" and t he lTurner et al.l (20051) 



" Heav y" component. We should note that both iTurner et al.l (120081 ) and IMehdipour et al. 
(I2OIOI ) used another component with an intermediate outflow velocity (~-1000 km s~^), 
namely "Zone 4" and component "C", respectively. Table 5 shows that the total column 
density can vary between observations and authors by nearly two orders of magnitude. The 
highest column densities are due to few spectral features e.g., Fe-K in the present work, 
or to the need to e xplain spectral curvature with photo-electric absorption without lines 
(iTurner et al.ll2008l ). The diversity in ionization parameter and velocity is mostly due to 
selective identification of spectral features. Our analysis shows that a broad range of both 
ionization and velocity are present. 



4.5. Possible Local Absorption 



The outflow velocity of component 3 is -2600±200 km s~^ matching the cosmological 
recession of -2650 km s~^, which raises the possibility that some of the absorption at this 
veloc ity is due to local absorption as common ly found along lines of sight to bright AGNs 



Nicastro et al 



2003 



Willianas et al.ll2005l ). A similar component was recently found in 



(e.g.. 

MCG -6-30-15 (iHolczer et al.ll2010l ). although there it was kinematically resolved (-1900±150 
km s~^ versus -2300 km s~^). The oxygen and nitrogen lines of component 3 are narrow 
and thus suspect of having a local origin. These lines require fturb = 100 km s~^ (FWHM = 
170 km s^^), which is less than fturb = 300 km s~^, which is used for higher ionization. This 
narro w width is consistent with local, ionized ISM UV absorption lines along this line of 
sight (IKraemer et al.ll2002l ). There are a few additional reasons to favor the local component 
scenario for O and N. First, they form at lower ionization parameters of log^ ~ — 2, 
while most of component 3 is primarily comprised of high ionization species. In Fe, we do 
not detect any M-Shell ions, and not even Fe"*"^^ in this component. Even though we can 
not conclusively determine whether the oxygen and nitrogen absorption comes from outflow 
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component 3, or has a local origin, we tend to favor the local origin scenario, because such 
a component is often observed. 

The (presumably) local component column densities are shown in Table O The low 
charge states are due to the Galactic disc and halo, and likely not associated with the higher 
ionization states that are due to the hot phase of the Galactic halo or the Local group. 
The oxygen ionic column densities are all in the range of several 10^^ cm~^, which implies 
a neutral hydrogen c olumn density of a few 10^*^ cm~^, using 4.9x10"^ for the O to H ratio 
( lAsplund et al.ll2009l ) and a fractional ionic abun dance of 0.5. This value i s consistent with 



the Galactic absorption of Nh = 3.23x10 cm (jPickey &: LockmadllQQOl ). Using the same 



calculations for N"*"^ gives a slightly higher neutral hydrogen column density of the order of 



40^1 cm-2. 



Note that most species of component 3 are likely not local since their ionization and 
column densities are too high. Ions such as Fe"'"^^""'"^^, K-shell S and Si, (see Table H]) are 
usually not observed in the local ionized ISM. Moreover, the high column density measured 
in this component of Nji^ ~ 10^'^ cm~^ is by far higher than typical local ISM columns. 



4.6. Variability 



The observation of 2006 caught NGC 3516 in a much higher state than the 2001 ob- 
servation, as can be seen in Fig. [2l So far, the analysis in this paper focused on the 2006 
observation. We now want to use the two flux states to study the differences between the 
two. Several e xplanations for the f lux and spectral variability of NGC 3516 can be found in 
the literature. iNetzer et al.l (120021 ) reported a slow, monotonic decay in flux between 1994 
and 2000, however BeppoSAX observations from 1996 and 1997 as well as the present data 
show a sharp transition b etween high and low flux over a period of 4 months (see Ta b le. [H) . 
Turner et al.l (l2005l . l2008l ) used a varying covering factor. However, iMehdipour et al.l (l2010l ) 
found that a varying covering factor did not fit the data, while, variable source continuum 
could. 

In order to compare the absorption in 2001 and in 2006 we scaled up the 2001 spectrum 
to the flux level of 2006 near specific lines. Since absorption depends exponentially on optical 
depth, this comparison directly shows changes in optical depth, independent of continuum 
flux. If the absorber did not change, the scaled-up trough should match that of the high- 
state. Because the continuum shape during the two states is different, this comparison is 
meaningful only locally. Indeed, we show this comparison around the most prominent lines. 
Results for Si and Mg K a lines, and for Fe"*"^^ and Fe"^^^ are shown in Fig. [TJ where the low 
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stat spectrum is multiplied by 4 in the upper left panel, by 4.5 in the upper right panel, by 
4.3 in the lower left panel, and by 5.2 in the lower right panel. The S/N in the low-state 
is worse, but it appears that component 1 (the slowest component) did not change much 
between the high and low states. On the other hand, some of component 2 and most of 
components 3 and 4 are absent in the 2001 low state. In the next section, we discuss a 
possible geometrical explanation for this result. 



4.7. Possible Geometry of Outflow 

We find that apart from a lower continuum level, the 2001 spectrum of NGC 3516 also 
lacks the faster absorption components. The appearance of high-ionization (log^ ~ 3.5 
erg s~^ cm) compo nents with columns of Nh ~ 10^^ cm~^ are reminiscent of the variable 



covering invoked by [Turner et al.l ( l2005l . 120081 ) to explain the varying continuum shape. We 



allude to three other possible explanations: 

A. Photo-ionization change - The faster components of the outflow are also more 
ionized. One possibility could be that the fast components recombined due to the reduced 
flux, and are thus not seen in 2001. However, one would still expect to observe the fast 
components in lower charge states. We could not detect the fast components in the 2001 
spectrum in any ion. Since the 2001 spectrum has a much lower signal to noise ratio, we can 
not unambiguously rule out this possibility. 

B. Fast components crossing line of sight - Another possibility could be that the 
fast components, while not in the line of sight in 2001, passed through our line of sight in 
2006, 5 years are plenty of time compared to the variability time scale of months that perhaps 
represents the size of the source, to make this scenario work. Such transverse velocities have 



been proposed for NGC 4151 by iKraemer et al.l (120061 ). 



C. Lighting up the disc The third possibility is that the fast components are present 
the whole time, but in 2001, when the flux is low, there was no light from accretion disk 
behind them for them to absorb. In 2006, the source flux is much brighter, whether intrinsic 
source output change or flux change due to variable covering factor by thick gas, which could 
result from a flare on the disc that illuminates the fast components from behind, and which 
subsequently they absorb along the line of sight. 

With the current spectra, one cannot rule out any of these scenarios. With better S/N 
spectra, one might be able to confirm (or rule out) changes in the (photo-) ionization state 
if the fast component is detected (or not detected) in low ionization species during the low 
flux state. If one would detect a change in the covering fraction in the narrow absorption 
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l ines, which we are n ot able to detect here, that would be evidence for transverse velocities 
( iKraemer et al.ll2006[ ). The possibility of flaring on the disc is most difficult to test, as the 
angular resolution required for detecting sources on sub-disc scales is currently prohibitive. 



5. CONCLUSIONS 

We have analyzed the kinematic and thermal structure of the ionized outflow in NGC 3516. 
We find absorption troughs in dozens of charge states that extend from zero to almost 5000 
km s~^. We model the outflow with four absorption systems. The flrst and second compo- 
nents are outflowing at -350 and -1500 km s~^, and span a considerable range of ionization 
from at least Fe"^^ to Fe^^^ [—0.5 < log^ < 3.5 (erg s~^ cm)]. The third and fourth compo- 
nents are outflowing at -2600 and -4000 km s~^, respectively and are highly ionized featuring 
only L-shell and K-shell iron and K-shell ions from lighter elements. Finally, a component 
of local absorption at z = (-2630 km s~^) is detected, this component could be part of the 
third component, however its low ionization and narrow time proflles imply it is more likely 
at z = 0. 

Using our AMD reconstruction method for all four components, we measured the dis- 
tribution of column density as a function of ^. We flnd a double-peaked distribution with a 
signiflcant minimum at 0.7 < log^ < 1.5 (erg s~^ cm) in components 1 and 2, which corre- 
sponds to temperatures of 4.5 < logT < 5 (K). This minimum was observed in other AGN 
outflows like MCG -6-30-15, NGC3783, IRAS 13349+2438, NGC7469, and it can be ascribed 
to thermal instability that appear to exist ubiquitously in photo-ionized Seyfert winds. The 
AMD of component 3 shows a continuous rise in column density toward higher ionization 
parameters. The AMD of component 4 is similar in some aspects to that of component 3, 
however it shows a minima at log T ~ 6 (K) mostly due to low Fe+^°-Fe^^^ column densities. 
The local absorption system could arise from either the ionized Galactic ISM, or from the 
local group. The fast components 3 and 4 are not present in the lower flux spectra of 2001. 

We thank Shai Kaspi for useful comments. This work was supported by a grant from 
the ISF. 



- 15 - 



REFERENCES 

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 
Badnell, N. R. 2006, ApJ, 651, L73 

Bar-Shalom, A., Klapisch, M., & Oreg, J. 2001, J. Quant. Spectr. Radiat. Transfer, 71, 169 
Behar E., Cottam J C, & Kahn S. M. 2001, ApJ, 548, 966 
Behar, E., Sako, M., & Kahn S. M. 2001, ApJ, 563, 497 
Behar, E., & Netzer, H. 2002, ApJ, 570, 165 

Behar, E., Rasmussen, A. P., Blustin, A. J., Sako, M., Kahn, S. M., Kaastra, J. S., 
Branduardi-Raymont, G., & Steenbrugge, K. C. 2003, ApJ, 598, 232 

Behar, E. 2009, ApJ, 703, 1346 

Blustin, A. J., Kriss, G. A., Holczer, T., Behar, E., & Kaastra, J. S. 2007, A&A, 466, 107 
Chelouche, D. 2008, larXiv:0812.362T] 

Costantini, E., Salvini, C., Comastri, A., Fruscione, A., Mathur, S., Nicastro, F., Stirpe, 
G. M., & Wilkes, B. 2001, X-ray Astronomy: Stellar Endpoints, AGN, and the Diffuse 
X-ray Background, 599, 590 

Crenshaw, D. M., Kraemer, S. B., Boggess, A., Maran, S. P., Mushotzky, R. F., & Wu, C.-C. 
1999, ApJ, 516, 750 

Crenshaw, D. M., Kraemer, S. B., George, I. M., 2003, ARA&A, 41, 117 

Decaux, V., Beiersdorfer, P., Osterheld A., Chen M., & Kahn, S. M. 1995, ApJ, 443. 464 

Detmers, R. G., Kaastra, J. S., Steenbrugge, K. C. et al. 2011, A&A, 534, A38 

Dickey, J. M., Lockman, F. J. 1990, ARA&A 28, 215 

Fukumura, K., Kazanas, D., Contopoulos, I., &Behar, E. 2010, ApJ, 715, 636 

Fukumura, K., Kazanas, D., Behar, E., & Contopoulos, I. 2011, Billetin of the American 
Astronomical Society, #327.17 

Gongalves, A. C, Goosmann, R. W., Mouchet, M., et al. 2010, Twlfth Inernational Solar 
Wind Conference, 1248, 439 



- 16 - 

Gu, M. F., Holczer, T., Behar, E., & Kahn, S. M. 2006, ApJ, 641, 1227 

Halpern, J., P., 1984, ApJ, 281, 90 

Haardt, F., et al. 2001, ApJS, 133, 187 

Holczer, T., Bchar, E., & Kaspi, S. 2005, ApJ, 632, 788 

Holczer, T., Behar, E., & Kaspi, S. 2007, ApJ, 663, 799 

Holczer, T., Behar, E., & Arav, N. 2010, ApJ, 

Kallman, T. R., & Bautista, M. 2001, ApJS, 133, 221 

Kaspi et al. 2001, ApJ, 554, 216 

Kaspi, S., Brandt, W. N., CoUinge, M. J., Elvis, M., & Reynolds, C. S. 2004, AJ, 127, 2631 
Keel, W. C. 1996, AJ, 111, 696 

Kolman, M., Halpern, J. P., Martin, C, Awaki, H., & Koyama, K. 1993, ApJ, 403, 592 

Kraemcr, S. B., Crenshaw, D. M., George, I. M., Netzer, H., Turner, T. J., & Gabel, J. R. 
2002, ApJ, 577, 98 

Kraemer, S. B., Crenshaw, D. M., Gabel, J. R., et al. 2006, ApJS, 167, 161 

Kriss, G. A., Espey, B. R., Krolik, J. H., Tsvetanov, Z., Zheng, W., & Davidsen, A. F. 1996, 
ApJ, 467, 622 

Kriss, G. A., et al. 1996, ApJ, 467, 629 

Krolik, J. H., McKee, C. F., Tarter, C. B. 1981, ApJ, 249, 422 

Krongold, Y., Nicastro, F., Brickhouse, N. S., Elvis, M., Liedahl, D. A., & Mathur, S. 2003, 
ApJ, 597, 832 

Krongold, Y., Nicastro, F., Elvis, M., et al. 2005, ApJ, 620, 165 

Krongold, Y., Nicastro, F., Elvis, M., et al. 2007, ApJ, 659, 1022 

Markowitz, A., et al. 2008, ApJ, 60, 277 

Mathur, S., Wilkes, B. J., & Aldcroft, T. 1997, ApJ, 478, 182 

Mehdipour, M., Branduardi-Raymont, G., & Page, M. J. 2010, A&A, 514, AlOO 



- 17- 

Nandra, K., Mushotzky, R. F., Yaqoob, T., George, I. M., & Turner, T. J. 1997, MNRAS, 
284, L7 

Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1999, ApJ, 523, 
L17 

Netzer, H., Chelouche, D., George, I. M., Turner, T. J., Crenshaw, D. M., Kraemer, S. B., 
& Nandra, K. 2002, ApJ, 571, 256 

Netzer, H., et al. 2003, ApJ, 599, 933 

Netzer, H. 2004, ApJ, 604, 551 

Nicastro, F., Zezas, A., Elvis, M., et al. 2003, Nature, 421, 719 

Piconcelli, E., Jimenez-Bailon, E., Guainazzi, M., et al. 2005, A&A, 432, 15 

Reynolds, C. S., 1997 MNRAS, 286, 513 

Sako, M. et al. 2001, A&A, 365, L168 

Sako, M. et al. 2003, ApJ, 596, 114 

Steenbrugge, K. C. et al. 2005, A&A, 432, 453 

Turner, T. J., et al. 2002, ApJ, 574, L123 

Turner, T. J., Kraemer, S. B., George, I. M., Reeves, J. N., & Bottorff, M. C. 2005, ApJ, 
618, 155 

Turner, T. J., Reeves, J. N., Kraemer, S. B., & Miller, L. 2008, A&A, 483, 161 
Williams, R. J., Mathur, S., Nicastro, F., et al. 2005, ApJ, 631, 856 



This preprint was prepared with the A AS macros v5.2. 



- 18 - 




Fig. 1. — Chandra HETGS spectra of NGC 3516 corrected for cosmological redshift {z = 
0.008836) and binned to 10 mA. Observation dates are listed in the legend. No absorption 
variability is observed on time-scales of days. 
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Rest Frame Wavelength (A) 



Fig. 2. — Chandra HETGS spectrum of NGC 3516 corrected for cosmological redshift {z = 
0.008836) and binned to 10 mA. The red and black hnes represent the average flux from the 
2006, and 2001 observations, respectively. 
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Rest Frame "Velocity (km s ') 



Fig. 3. — Chandra HETGS 2006 average spectrum of NGC 3516 in velocity space, corrected 
for cosmological redshift {z = 0.008836) around Si+^\ Si+^'^, Mg+^\ Mg+i°, Fe+^s, Fe+^\ 
Fe"*"^^ and Fe"*"^^ lines. Four kinematic components can be discerned. Other lines are marked 
in green as are two unidentified troughs next to Mg"*"^^ (see also Figure H] at 8.2 A). For Fe"*"^^ 
only the two slow components are present, which is an indication of the high ionization of 
the fast components. 
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Fig. 5. — AMD of component 1 (solid black) and component 2 (dashed blue) in the outflow of 
NGC 3516, obt ained exclusively from Fe absorption and scaled by the solar Fe/H abundance 
of 3.16x10"^ (lAsplund et al.l 120091 ). The corresponding temperature scale, obtained from 
the XSTAR computation is shown at the top of the figure. The accumulative column density 
up to ^ is plotted in the lower panel, yielding a total of Nh = (1-8 ± 0.5) x 10^^ cm~^ for 
component 1 and Nh = (2.5 ± 0.3) x 10^^ cm~^ for component 2. 
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Fig. 6. — AMD of component 3 (black) and component 4 (dashed blue) in the outflow of 
NGC 3516 obt ained exclusively from Fe absorption and scaled by the solar Fe/H abundance 
of 3.16x10"^ (lAsplund et al.l 120091 ) . The corresponding temperature scale, obtained from 
the XSTAR computation is shown at the top of the figure. The accumulative column density 
up to ^ is plotted in the lower panel, yielding a total of Nh = (6.9 ± 4.3) x 
component 3 and Nh = (5.4 ± 1.2) x 10^^ cm~^ for component 4. 
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Fig. 7. — Chandra HETGS spectrum of NGC 3516 corrected for cosmological redshift {z 
= 0.008836) around K-Shell Si (upper left panel), K-Shell Mg (upper right panel), Fe"*"^^ 
(lower left panel) and Fe"*"^^ (lower right panel) lines. The black and red spectra represent 
the combined observations of 2001, and of 2006 respectively. The blue data represent the 
2001 spectrum scaled up to match the 2006 continuum. The four kinematic components are 
labelled. The spectral lines in the 2001 lack absorption in the fast component that is present 
in 2006. The lower panels are the ratios between the 2006 and the 2001 scaled spectra. 
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Table 1. Historical observations of NGC 3516 



Observatory 


Year 


Duration Flux level at 1 keV 


References ^ 






(ks) (10-3 ph cm-2 s-^ keV-i ) 





Ginga 


October 1989 


20 


6±1 


1 


ROSAT 


October 1992 


13 


17±3 


2 


ASCA 


April 1994 


28 


^13±2 


3 


ASCA 


March 1995 


33 


^7±1 


4 


BeppoSAX 


November 1996 


16 


1±0.5 


5 


BeppoSAX 


March 1997 


16 


7±2 


5 


ASCA 


April 1998 


360 


^3±0.5 


6 


Chandra 


October 2000 


47 


^0.6±0.2 


6 


Chandra 


April 2001 


119 


1.5±0.2 


7,11 


XMM-Newton 


April 2001 


49 


1.7±0.2 


7 


Chandra 


November 2001 


88 


0.9±0.2 


7,11 


XMM-Newton 


November 2001 


49 


1±0.2 


7 


Suzaku 


October 2005 


135 


0.6±0.2 


8 


Chandra 


October 2006 


200 


5±0.5 


9,11 


XMM-Newton 


October 2006 


155 


5±0.5 


9,10,11 



'^l. Kolman et al. ( 


199i. 2. Mathur et al.l ( 


1997). 3. Reynolds (1997). 4. Kriss et al. 


Costantini et al. 


(2000 


). 6. Netzer et al. ( 


2002). 7. Turner et al. ( 


2005). 8. Markowitz et al. 


Turner et al. ( 


20081. 10. 


Mehdipour et al. 


(2010tl 


11. present work. 



iNetzer et al.l (j2002l) finding of long-term decline over years is fortituous. 
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Table 2. Chandra observations of NGC 3516 used in this work 



Obs. ID 


Start 


Detector 


Gratings 


Exposure 


Counts in HEG 


Counts in MEG 




Date 






(s) 


orders ± 1 


orders ± 1 


2431 


2001 April 9 


ACIS-S 


HETG 


35568 


2383 


4143 


2080 


2001 April 10 


ACIS-S 


HETG 


73332 


9019 


16057 


2482 


2001 November 11 


ACIS-S 


HETG 


88002 


6393 


11181 


8452 


2006 October 9 


ACIS-S 


HETG 


19831 


6245 


12024 


7282 


2006 October 10 


ACIS-S 


HETG 


41410 


9176 


16906 


8451 


2006 October 11 


ACIS-S 


HETG 


47360 


17482 


34473 


8450 


2006 October 12 


ACIS-S 


HETG 


38505 


14773 


29920 


7281 


2006 October 14 


ACIS-S 


HETG 


42443 


9572 


16959 



Table 3. Narrow Emission Lines ^ 



Line 


•^Rcst 


\ b 

^Observed 


Flux 




(A) 


(A) 


(10~^ photons s~^ cm~^) 


Fe+o - Fe+9 Ka 
Fe+io - Fe+16 Ka 


1.94 
1.93 - 1.94 ^ 


1.936 ± 0.01 


3.3 ± 0.6 


Ne"*"^ forbidden 


13.698 


13.69 ± 0.01 


0.5 ± 0.1 


0+^ forbidden 


22.097 


22.093 ± 0.01 


6 ± 1 



^FWHM = 235 km s^^ applied uniformly to oxygen and neon emission lines. FWHM 
3500 km was applied to iron Ka emission line 

'^in the AGN rest frame. 



Decaux et al.l (|l995l ). 
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Table 4. Current best-fit column densities for ions detected in the 2006 HETGS spectrum 

of NGC 3516. 



Ion Column Density 



'ion 

(IQis cm-2) 



Column Density 

Nion 

(IQie cm-2) 



Ion Column Density Column Density 



N- 

(10^6 cm-2) 



(IQis cm-2) 



Compl Comp2 Comp3 Comp4 



N+6 



O+o 
0+1 

0+2 
0+3 
0+4 
0+5 
0+6 

0+7 



Ne+3 
Ne+^ 
Ne+5 
Ne+6 
Ne+7 
Ne+8 



Mg+4 

Mg+5 

Mg+<5 

Mg+7 

Mg+8 

Mg+9 

Mg+io 

Mg+" 



Si+5 
Si+6 
Si+7 

Si+« 

Si+9 
Si+10 

Si+" 
Si+12 

Si+13 



50 



3.01?:? 



+3 

11 
+ 1.1 
5.0 
+ 1.4 
-3.7 
+90 
1.0 



11 

5.0 
4.0 
10 

18+300 



-4.0 
+2.0 
6.0 

2.011^ 

+2.7 



4.0 
6.0 



2.0 



-2.0 
o n+0.7 
"^■'J-S.O 
+0.4 
-2.0 
+80 
2.0 



-1.0 

-,+1.7 

-2.5 
+ 1.3 
-0.9 

4 0+''-^ 

^■^-0.8 



2.0 
20 

100t?o° 



1.0^ 



4.0Z 
1.0 



4.0 



+ 18 



-0.4 
-:',() 
-0. 1 



*-^-2.0 
0.5t3-2 

3.0 

2.0 
2.0 
6.0 



u-'J-o.s 

■■■•3-0.4 

1 t;+3.0 
J^-3-0.2 

2 0+2-^ 
2 0+1-^ 
4 0+^-* 



0.5 
+0.8 
-2.5 
+6.0 
0.4 
+ 1.1 
1.0 
+ 25 
-0.6 
-.-,() 

l.;i 



5.5 



+ 1.9 
-0.6 



'^•'^-0.8 
o n+3.2 
''•'-'-0.8 

■^•^-0.3 
o n+2.8 
'^■'-'-0.3 

2 0+1-2 

^•'-'-0.4 

2 6+° *5 

^•"-0.5 
9 c+O.S 

^•°-0.7 
4 7+4.3 
^•'-0.4 



^•"-0.6 
1 c+3.5 
■'■■''-0.2 

4 0+2" 

^■"-0.4 
^•"-0.4 



^+2.0 

'-'-0.2 

-0.7 



4 

^•^-0.8 

o n+2.3 
"5-"-0.6 

^•"-0.3 

^•^-0.3 

+ 1.8 
-0.7 
+7.7 
1.5 



g+14 

S+15 



Ar+16 
Ar+17 



Ca+18 

Ca+19 



18 
0.6 



Q+3.2 
^-0.1 
c+2.8 
''-O.S 



Fe+1 
Fe+2 
Fe+3 
Fe+4 
Fe+5 
Fe+6 
Fe+7 
Fe+8 
Fe+9 
Fe+io 

Fe+ii 
Fe+12 
Fe+13 
Fe+14 
Fe+15 
Fe+16 
Fe+17 
Fe+18 
Fe+19 
Fe+20 
Fe+21 
Fe+22 
Fe+23 
Fe+24 
Fe+25 



7.0 
15 



-+2.5 



14+6.: 



i-5t°o:^ 



Compl 

-1-^2.3 
-0.9 
+4.1 
-1.8 



2.o: 

6.0 



Comp2 

2.1 



3.0] 
5.0 



3.0 
3.0 



1.7 
+4.3 
2.1 



+1.9 
-3.6 



Comp3 Comp4 

rFTH — 

-2.0 



3.5 



8.0 



Q+6.0 
^-1.8 



2.0 



+6.6 

-2.0 



U.Z_Q 2 

0.2l°i 
0.2l°i 
0.2t°i 

^•'^-0.2 

5+21 

1 0+1° 

1 0+^1 
■■■•^-0.2 

5+1-7 

U.O_Q 4 

-I 7+0.8 
4-'-0.7 

o.slgj 

0.21°;^ 

2 4+0.5 

^•^-0.5 
9 c+0.2 

5 

o p + 1.0 
0.0_Q 4 

9 9+2.2 
^•^-0.2 

3.0+3-° 



4 0+^^-° 



5.01°; 



-0.3 

3-oli^:^ 
o.3l«:^ 

o.ii^"i 
o.it^^i 



0.2: 



-0.2 



-,+2.6 
-'-2.9 



0.2 
2.0 
2.0 



o.2i°:i 
o.2i°:3 

2+°-'i 

2+"-3 

+0.3 
0.2 
+ 1.7 
0.3 
+2.7 
0.3 

^■^-0.7 

A 0+1-7 

1 4+1-2 

9 f.+0.4 

1 o+°-3 

■■■•^-0.4 

o.5l°:2 
o.5tJ:? 

9 c + 1.6 
9 c+0.6 

e.olJ:^ 

3 q+2.4 

•5 c+8.8 
o.0_o.4 

o r+4.5 
o.O_o 4 

o Q+7.9 

-12 
-0.4 
+ 70 
-12 

10+7° 



4.0: 

30 



1+7.0 
-1.5 



3.01^7 
4.4 

4.2 



+0.8 
0.7 
1.1 
-0.5 
+9.6 
0.3 



2.8 

2.71^:^3 

4.5 
0.3 
+4.5 
-0.8 
+66 
14 
+ 150 
20 



3.0 

8.0 
22 
60 



1 5+1-6 

^■■''-0.2 

1.3: 

2.0 



1.0 



0.8 
0.4 
0.6 
-0.2 
+2.4 
-0.2 



0.71^:1 

o.5tJ:^ 



hO.6 
-2.0 

h25 

130+1° 



2.0: 

20^ 
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Table 5. Physical parameters for absorption components in order of outflow velocity: 

comparison. 



Ref. 


Observatory 


Component 


Outflow 


Column 


Total Column 


Ionization 








Velocity 


Density 


Density 


Parameter 








(km s-^) 


(10^1 cm-2) 


(10^1 cm-2) 


lofff fere cm.) 


1 


ASCA 


1 

2 


<-120 
<-120 


14.1±3.9 
6.9±0.6 


21±4 


1.66±0.31'' 
0.32±0.11'' 


2 


ASCA 








10±1 


1.44±0.03 


3 


ROSAT 




-500 


7±1 


7±1 


0.90-1.11 ^ 


A(i 
4a 


BeppoSAX 


Warm 
Hot 


-500 
-500 


10±0.4 
158±73 


168±73 


0.73±0.10'^ 
2.36±0.10'' 


c;e 



AdLA 






Q 


Q 


b 
-Z.D 


D 


XMM-Newton 
Chandra 


UV 
High 


-200 

-ilOO 


6±2 

1 a 
10 


OTO_LOO 
Z/Z±ZO 


-0.5 

c 
z.o 




Heavy 


1 1 f\f\ 
-iiOO 


zoO±zo 




0.0 


7 


Suzaku 


Primary 
High Ion. 

Zone 1 


-ilOO 


55±2 

/in+46 
40_3i 

2 4+0.3 
^■^-0.2 


95±45 


0.3±0.1 

q 7+0.3 
-J-' -0.7 

9 /iq+0.58 

z.'±0_o,o3 


8 

(J 


XMM-Newton 


Zone 2 




0.5±0.1 


t: w 1 _l — L J. W 


0.25 


Chandra 


Zone 4 


-1000 


262+^^ 


4 Q-i+1.19 






Zone 3 


-1600 


202l«^ 




2.19±0.07 






A 


-100 


4 




0.9 


9 


XMM-Newton 


C 


-900 


10 


34 


3.0 






B 


-1500 


20 




2.4 






Comp 1 


-350±100 


2.8±0.6 




-0.5 - 0.7 






-350±100 


15.3±5.0 




1.5-3.5 


10 


Chandra 


Comp 2 


-1500d=150 

-1500±150 


5.9±0.6 
18.8±3.5 


166±45 


-0.5 - 0.7 
1.5 - 3.5 






Comp 3 


-2600±200 


69±43 




2-4 






Comp 4 


-4000±400 


54±12 




2-4 
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' ^ijKriss et a l.' ^199^, 2 lRevnoldsl JiQQtI). siMathur et alj JiQQtI) . 4lCos"tantini et alj (|2000|). 5. iNetzer et all |2002l ). 
6. [Turner et al- (,2005,), 7. iMarkowitz et all |2008h . 8. iTurner et all (|2008l ). 9. iMehdipour et all (|20ld) 10. present 



work. 

'^Ionization parameter used is U. 
'^Ionization paramer used is Vqx- 
'^Reference to data set 97BF therein. 



"Reference to data set ASCA98 therein. 
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Table 6. Oxygen and nitrogen ionic column densities at z — 0. 



Charge 


A^ion local 




State 


(lO^*' cm-2) 


(A) 


N+6 


20+f 


19.825,20.911 


neutral 


4.0+f 


23.523 






23.347 


0+2 




23.071 


0+3 


3.01^:^ 


22.741 


0+4 


-'-•'-'-1.5 


22.374 


0+5 


-'-•'-'-1.0 


22.019 


0+6 


^-0.5 


21.602 


0+^ 


lOlf 


18.969 



